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A generalization of the random geometric graph (RGG) model is proposed by considering 
a set of points uniformly and independently distributed on a rectangle of unit area instead of 
on a unit square [0,1]^ . The topological properties of the random rectangular graphs (RRGs) 
generated by this model are then studied as a function of the rectangle sides lengths a and 
b = 1/a, and the radius r used to connect the nodes. When a = 1 we recover the RGG, and 
when a —?■ oo the very elongated rectangle generated resembles a one-dimensional RGG. We 
obtain here analytical expressions for the average degree, degree distribution, connectivity, 
average path length and clustering coefficient for RRG. These results provide evidence that 
show that most of these properties depend on the connection radius and the side length 
of the rectangle, usually in a monotonic way. The clustering coefficient, however, increases 
when the square is transformed into a slightly elongated rectangle, and after this maximum 
it decays with the increase of the elongation of the rectangle. We support all our findings 
by computational simulations that show the goodness of the theoretical models proposed for 
RRGs. 

TAGS: 89.75.-k; 02.10.Ox 


INTRODUCTION 


The use of graphs for representing physical systems is becoming ubiquitous in many areas 
of theoretical and applied physics [ij. We can mention the use of graphs in statistical me¬ 
chanics and condensed matter physics, for solving Feynman integrals as well as in the study 
of quantum phenomena Q a. Mo. .ceot,. the u. of .as been ves. b.adeaed 

by their application in the analysis of complex systems 3|-l5|. In this case, those graphs 
receive the name of complex networks, due to the fact that they represent the skeleton 
of complex interconnected systems. In this case, networks are used to study a variety of 
physical scenarios, ranging from social and infrastructural, to biological and ecological ones. 
Here, we will use the terms graphs and networks interchangeably. When graphs are used to 
represent real-world physical systems it is necessary to have at our disposal some null model 
that allows us to evaluate which properties of the system have arisen from their connectivity 
pattern. In this sense, the common election is the use of random graphs. These are graphs 
with the same number of nodes and edges as the one under study, but in which the connec¬ 
tion between the nodes is made randomly and independently j^. There are several of these 
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random models o^reat usability in current network theory, such as the Erdos-Renyi j^, the 
Barabasi-Albert 9j or the Watts-Strogatz model to mention just three. 

In many real-world scenarios the networks e merg e under certain geometrical constraints. 


This is the case of the so-called spatial networks ll|, which include infrastructural networks 
such as road networks, airport transportation networks, etc., ll| and certain biological 
networks such as brain networks or the networks representing the proximity of cells in a 
biological tissue (see j^). The list also includes the networks of patches and corridors in 


a landscape 


121 , 1 


le networks of galleries in animal nests 


fractures in rocks 


Q,Qi, 


and the networks of 


ISl - among others. The classical election of a random graph used to 
represent these systems are the so-called random geometric graphs jlG, 17|. Here the term 
random geometric graph (RGG) is reserved for the case in which the nodes of the graph are 
distributed randomly and independently in a unit square and two nodes are connected if 
they are inside a disk of a given radius centered at one of the nodes. Other graphs in which 
the edges are constructed by using different geometric rules will be named here generically 
as random proximity graphs. 

RGGs have found important applications in the area of wireless communication devices 
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20| . such as mobile phones, wireless computing systems, wireless sensor networks, etc. 


This was indeed the hrst application in mind when Gilbert proposed the very hrst RGG 


model 


2l| . RGGs have also found applications in areas such as modeling of epidemic 


spreading in spatial populations, which may include cases such the spreading of worms in 
a computer network, viruses in a human population, or rumors in a social network 
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26| . RGGs have been used to describe how cities have been evolving under the geometric 
constraints imposed by their geographic locations j^|. For a wider perspective on the 
applications of spatial graphs the reader is referred to the review {h|. 

In all the previously mentioned real-world scenarios, the shape of the location in which 
the nodes of the graph are distributed may play a fundamental role in the topological and 
dynamical properties of the resulting graphs. That is, it is intuitive to think that the con¬ 
nectivity, distance, clustering and other fundamental topological properties of the graphs 
are affected if we, for instance, elongate the unit square in which the points are distributed. 
Here, we develop a new model that generalizes the RGG by allowing the embedding of the 
nodes in a unit rectangle instead of a unit square. Our main goal is to investigate how 
the elongation of a unit square influences the topological properties of the graphs generated 
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by the model. These generalized graphs will be named here the random rectangular graphs 
(RRGs). In this work we concentrate on the influence of the length of the rectangle on the 
topological properties of the graphs emerging on them, such as their average degree, connec¬ 
tivity, degree distribution, average path length and clustering coefficient. In particular, we 
And analytical expressions and bounds for all of them and provide computational evidence 
of the goodness of these approaches for relatively large RRGs. 


II. DEFINITION OF THE MODEL 


The RGG is defined by distributing uniformly and independently n points in the unit d- 
dimensional cube [0,1]*^ |16| . Then, two points are connected by an edge if their (Euclidean) 
distance is at most r, which is a given fixed number known as the connection radius. 


Let us now define a unit hyperrectangle as the Gartesian product [oi, hi] x [ 02 , 62 ] x • • • x 
[orf, where a*, 6* G M, a* < hi, and 1 < i < d. Hereafter we will restrict ourselves to 
the 2 -dimensional case, which corresponds to a rectangle of unit area, which we will call 
the unit rectangle. Now, the RRG is defined by distributing uniformly and independently 
n points in the unit rectangle [a, h] and then connecting two points by an edge if their 
(Euclidean) distance is at most r. It is evident that the only change we have introduced 
here is to consider a rectangle of unit area instead of the analogous square. The rest of the 
construction process remains the same as for the RGG. This means that RRG —)■ RGG as 
(a/ 6 ) ^ 1. In this sense we can say that the RRG is a generalization of the RGG. In Fig. [T] 
we illustrate an RGG and an RRG constructed with the same number of nodes and edges. 


An interesting question is what happens at the other extreme, when a —)■ cxo. In this 
case we have that 6 —)■ 0 , which means that the n points are uniformly and independently 
distributed on the straight line. Let us now consider a disk of radius r > 0 centered at each of 
these points and let us connect every point to the other points which he inside its disk. Thus, 
the resulting graph resembles a one-dimensional RGG, that is a graph created by placing 
the n points uniformly and independently on the interval [ 0 , 1 ] and then connecting pairs 


of nodes if they are at a ( 
radius r (see for instance 


jUC 


28 


idean) distance smaller than or equal to a certain connection 


2!l). 
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(a) 



(b) 

Figure 1. Illustration of two random rectangular graphs with a = 1 (top), which corresponds to 
a random geometric graph on a unit square and with a = 2 {h = 0.5) (bottom). Both graphs are 
built with 500 nodes and 1750 edges. 

III. TOPOLOGICAL PROPERTIES OF RRGS 

A. Average degree 

We start the study of the topological properties of RRGs by considering an analytical 
expression for the average degree k. We remind that the degree of a node is the number 
of edges connected to it. The average degree is a property not only interesting by itself 
but as we will see in the next sections of this paper it is fundamental to understand other 
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topological parameters of RRGs. 

To start with, let us consider that for a given node, there are n — 1 nodes distributed in 
the rest of the rectangle. Dehne Ap to be the area within the radius r of a point p which 
lies within the rectangle. Since the nodes are uniformly and independently distributed, the 
expected node degree of a node Vi is E(/cj) = (n — \)Ai/{ah), where Ai is taken for the point 
where node Vi is located. This is because dividing the nodes between the area within distance 
r and the rest of the rectangle gives rise to the Binomial distribution Bin{n — Ai/{ah)). 
Averaging this over all possible node locations (i.e., the points in the rectangle) gives 

_ fp((^ - -1) Ip A . 

~ ah ~ (afe)2 ■ ^ ^ 

Let /(a, h, r) be the area within radius r of a point which lies in the rectangle, integrated 
over all points, i.e., f{a,h,r) = f^Ap. Based on preliminary computational results (not 
shown) obtained for the average degree we consider here the following three regions: 0 < 
r<h, h<r<a and a < r < \/a^ + 6^, recalling that a > h. We call these cases 1, 2 and 
3, respectively. Thus, the function f{a,h,r) takes different forms fi for each case i. This 
means that we can write 


with 


(n — 

(a6)^ 



( 2 ) 


h=< 


fi 0 < r < 6, 

/2 h<r <a, 

/s a <r < y/aF+lA, 


( 3 ) 


and our task is now to hnd the analytical expressions for fi. 

We consider the rectangle in Fig. [21 which shows 3 quarter circles of different radii (each 
corresponding to one of the three cases) as they intersect the interior of the rectangle. We 
consider only quarter circles instead of circles, then quadruple the result at the end. For 
each of these quarter circles, we divide them into vertical rectangular strips of width Ax 
which will approximate the areas of the intersection between the quarter circles and the 
full rectangle; this will become exact in the limit as Ax 0. We now consider several 
possibilities for these strips. 
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Figure 2. Illustration of three different quarter circles in the rectangle corresponding to 0 < r < 6 
(solid line), b < r < a (broken line) and a < r < y/a? + h‘^ (dotted line). The direction used for 
displacing the circles is represented as bottom-top (B-T) with an arrow in the graphic. 


First, the strips may approximate an area which is not rectangnlar, which occurs when 
the height of the strip is smaller than the height b of the rectangle. For a strip of distance 
X from the left of the rectangle, this corresponds to 0 < x < r for the smaller quarter circle 
(case 1), y/r^ — < x < r for the medium quarter circle (case 2), and y/r^ — b'^ < x < a 

for the largest quarter circle (case 3). Setting p = min{a,r), q = min{b,r), we have 

y/r^ — < X < p. 

Since we need to integrate these areas over all possible quarter circles, given a fixed 
radius, we wish to know how we can translate the quarter circle in the hgure and preserve 
a particular strip. That is, for a particular strip of distance x from the left of the rectangle, 
we may hnd a corresponding strip on the other quarter circles of the same radius. Since we 
have a rectangular strip of width Ax, height y/r^ — x^, and distance x from the center of 
the (quarter) circle, we may hnd this strip in any of the (a — x — Ax) positions horizontally, 
and (b— y/r"^ — x^) vertically. Thus, we can use integration to hnd the total area of all these 
strips by multiplying (a — x — Ax) (6 — y/r^ — x^) by the area of the strip y/r^ — x^ Ax, and 
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taking the limit to obtain 


h 



— dx 



— x^ — 


(r^ — x^)) dx. 


( 4 ) 


Secondly, we note that if these strips are translated far enongh in the bottom-top (B-T) 
direction, they become trnncated by the top of the rectangle. For a particular truncated 
strip, we may still find a corresponding strip in any of the (a — x — Ax) positions horizontally, 
and the truncated height t of a strip may be any value between 0 and the full height of the 
strip. Thus, integrating gives 


I 2 



r*y/r 


X 


t dt dx 



x^) dx. 


( 5 ) 


Alternatively, we may have \/r^ — 6^ > 6, in which case the rectangular strip is exact 
and of height b. In this case, the only contribution is from the truncated strips. We note 
that this applies for 0 < x < by a similar argument as before, and we integrate as 

follows 


rb 

Is = / {a — x) tdtdx 

Jo Jo 




.^1 


-(a — x)b‘^ dx. 


( 6 ) 


Thus, we have the expression for / as four times the sum of the above integrals 


/ = 4(Ji + I 2 + h) 


-Vr^-q 


2{a — x)b‘^dx 




_(a — x){AbVr‘^ — — 2(r^ — x‘^))dx, 


( 7 ) 


which can be written as 

















0 < r < 6 TTr'^ab — |(a + b)r^ + 

b <r < a —|ar^ — + a(|r^ + |6^)\/r2 — 

+2r‘^ab arcsin(^), 

a <r < ^/oF+lP —r‘^{o? + 6^) + |(a^ + 6^) — 

+6(|r^ + |a^)\/r2 — + a(|r^ + |6^)\/r2 — 6^ 

—2a6r^(arccos(^) — arcsin(^)). 


(8) 


Now we evaluate computationally how good the expression for the average degree of the 
RRGs is. In Fig. |3t^a) we plot the values of the average degree observed for RRGs with 
three different values of the rectangle side length. These observed values (represented by 
solid squares, circles and triangles) are the average of 100 random realizations of RRGs with 
n = 1, 500 nodes. The solid lines represent the expected values according to the expressions 
(I8|) . The Pearson correlation coefficients for the linear regression between the observed and 
expected values are larger than 0.9999 in the three cases. We enlarge the region of small 
radii for the case a = 30 (see Fig. [3] (b)) where it can be seen that it is a perfect ht also for 
this region with Pearson correlation coefficient as good as for the general case. 


B. Connectivity 


In this section we are interested in determining the connection radius for specihc values 
of n, i.e., r (u), which guarantees that the RRG F (n, r (u)) is asymptotically connected with 


3l| proved that if is the maximum length of an 


probability one. For the RGG, Penrose 
edge in the graph, then the probability that nnM^ — logn < a for a given a G M is 


lim P \mrM^ — logn < a\ = exp (—exp (—a)). (9) 

n^oo *- -* 

In 2D Mn = r , such that we have for the RGG 


lim P Irnir'^ — logn < a] = exp (— exp (—a)). (10) 

This means that for a —?> +cxd the RGG is almost surely connected when n —?■ oo, and 
almost surely disconnected when a —)■ — cxo. 
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Figure 3. (color online) (a) Illustration of the fit between the observed (black squares (o = 1), red 
circles (a = 3), and green triangles (a = 30)) and expected (solid line) values of the average degree 
for RRGs with different side lengths of the rectangle, (b) Wider range of radii for a = 30 (zooming 
for the small radii in the inset). 


The term is just the average degree k in RGG (when no boundary effects are 

considered). Thus we can write 


lim P \k — \ogn < «1 = exp (—exp (—a)). (11) 

In the previous section we have obtained an analytic expression for k in the RRG. Thus, 
we can replace this in (ITT]) and obtain an analogous expression for the RRG. That is, for a 
RRG we have that 


lim P 

n—)-oo 


{n - l)/i 
{ahf 


— log n < a 


exp (— exp (—a)), 


( 12 ) 


where fi is given by flHl) . 

Because the parameter a is unknown and it depends on the specihc RRG considered, we 
can obtain a lower bound for exp (—exp (—a)) using flT^ . That is. 
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exp 


< exp (— exp (—a)). 


(13) 



({n- l)fi 

V (aby 



Now we consider the computational evaluation of the connectivity of RRGs with 1,500 
nodes as a function of the connection radius and the rectangle side length. We start by 
analyzing the goodness of the upper bound that we have found for the probability of a RRG 
to be connected (see (IT^ h In Fig. |l](a) we illustrate the plot of the probability of being 
connected versus the connection radius for graphs embedded into rectangles of side length 
a = 1 (black circles) and a = 10 (red circles). The solid circles represent the average values of 
100 random realizations for these graphs. The values corresponding to the upper bound are 
plotted as broken lines. As can be seen, both the observed and the upper bound follow the 
same distribution and the upper bound is relatively close to the average observed values. In 
Fig. |4](b) we illustrate the change in the probability that the RRG is connected as a function 
of the connection radius, i.e., the minimum radius needed to make the network connected, 
for three different values of a. As the square is elongated the critical radius increases with 
the value of a. For instance, for a = 1 the critical radius is about 0.2, and for a = 10 it is 
about 0.9. Then, for a = 1, a > 185.3 for the RRG to be connected. This value increases 
up to a > 750.8 for a = 5 and to a > 3813.8 for a = 10. The main reason for this increase 
in the critical radius is that as we elongate the rectangle the points have to cover a longer 
region and consequently their separation increases. As a consequence, we need to increase 
the connection radius in order to guarantee the connectivity of the network. In other words, 
increasing the value of a necessarily implies having to increase the connection radius to make 
the network connected. The global increase of the critical radius with the side length of the 
rectangle is given in Fig. H] (c). 


C. Degree distribution 


In the RRG the n nodes are distributed uniformly and independently on the unit rect¬ 
angle. Then, the degree distribution can be easily estimated by considering the probability 
density function of having a node i of degree ki given that there are n — 1 other nodes 
uniformly distributed in the unit rectangle (see for instance 32j). This gives rise to the 


binomial distribution of the node degrees, which when (n — 1) ~ n is given by 
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(a) (b) (c) 

Figure 4. Dependence of the probability of being connected with the connection radius in RRGs. 
a) Illustration of the goodness of the bound (1131) for two RRGs with a = 1 (black curves), and 
a = 5 (red curves). The observed values for RRGs with 1,500 nodes averaged over 100 random 
realizations are illustrated as solid circles and the values obtained by the theoretical bound are 
given by broken lines, (b) Observed values of the probability of being connected for RRGs with 
1,500 nodes averaged over 100 random realizations and three different rectangle side lengths a = 1 
(black squares), a = 5 (red circles), and a = 10 (green triangles), (c) Variation of the radius for 
which the RRG is connected (critical radius) with the rectangle side length for RRGs with 1,500 
nodes and 11,250 edges. 


p(k) 


^ j / (1 - p) 


(14) 


When n — oo and p is sufficiently small, this binomial distribution approaches very well 
a Poisson distribution of the form 


p (k) ~ 


k^ exp (—k) 
k\ 


( 16 ) 


As we have previously obtained an analytic expression for the average degree k we can 
easily compute the degree distribution for RRGs. We select RRGs with 5,000 nodes and 
radius of connection equal to 0.025. Then, we obtain the degree distribution for different 
values of the rectangle side length and take the average of 100 random realizations. In Fig. [5] 
we also plot the expected distribution using the equation (IT^ in which we have plugged the 
values of the expected average degree obtained previously. As can be seen, independently 
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0.35 


0.3 



k 

Figure 5. Degree distribution of RRGs with n = 5,000, connection radius r = 0.025, and rectangle 
side lengths a = 1 (triangles), a = 50 (circles), a = 100 (squares). The solid figures (triangles, 
circles and squares) correspond to the average of 100 random realizations for the given network. 
The solid lines correspond to the shape of the Poisson distribution (fT5ll with the corresponding 
average degree obtained from eq. ([8]). 

of the side length of the rectangle, the RRG displays Poisson degree distributions. That is, 
the elongation of the rectangle does not affect the shape of the degree distribution of the 
nodes. 


D. Average shortest path distance 

Let r = (V, E) be a simple connected graph. A path of length /c in P is a set of nodes 
... ,ik,ik+i such that for all 1 < / < A;, G E with no repeated nodes. The 

shortest-path or geodesic distance between two nodes u, n G R is dehned as the length of the 
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shortest path connecting these nodes. We will write d{u, v) to denote the distance between 
u and V. Here we will define, as is usual in network theory, the average path length to be 
the following quantity: 


= ( 16 ) 

^ ^ u<v 

The diameter of the graph is defined as the maximum of all the shortest path lengths in 
the graph, i.e., D = max^ ^d {u,v). We consider here an upper bound for the average path 
length. Then, we start by considering that the n points are distributed homogeneously in 
the rectangle of sides a and b = a~^ in such a way that the points are equally spaced in the 
rectangle and are separated by a Euclidean distance r. In this case, the point which is the 
farthest from the rest of the n — 1 points in the rectangle is one at any of the four corners of 
the rectangle. Let us designate this point as i. If the average path length of this node is (/j), 
it is straightforward to realize that {/) < (/,), and we will have the desired upper bound. It 
is easy to see that the largest chain of connected points is realized along the main diagonal 
of the rectangle. Thus, the longest distances involving the node i are those connecting it 
with the other nodes j along the main diagonal of the rectangle. Let us designate the length 

Q 

of the main diagonal is c. Then, there are D = - connected nodes in this line. Let us now 

r 

consider (h) based only on those j nodes along the main diagonal of the rectangle (notice 
that this is an upper bound for (/*)). That is 


Consequently, 


j£diag 


1 + 2 + --- + D 


D{D + l) 
2 


(y < = ELI 

^ ^ - 2D 2 

Because, (/) < (/,), we have 


(0 < ih) < 


D + 1 
2 


It is easy to note that D < — -. That is, the length of the diagonal is smaller than 

r ar 

the sum of the length of the two sides of the rectangle (in terms of number of nodes), which 

are — and - = —. Then, we have 
r r ar 
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a 1 

D + 1 ~ - 

(/) < ik) < 


+ ar + 1 


(17) 


2 ~ 2 2ar 

An important consequence of this bound is that it allows us to find the asymptotic 
behavior of the average path length as the rectangle becomes more elongated. That is, 
if we fix the connection radius we can see that linia^oo (0 = other words, as the 

rectangle becomes more elongated the RRG becomes a large-world with a very large average 
shortest path distance. As a —)■ oo, the connection radius r has to be increased to keep the 
connectivity of the RRG. Thus, the infinite growth of the path length is not observed in 
practice due to the connectivity constraints imposed by the connection radius. Our further 
calculations show that the average path length is of the same order of magnitude as the 
right-hand side of eq. (fT7|) . which means that we can write it as: 


(0 


-I- ar -I- 1 
2ar 


( 18 ) 


We now study computationally RRGs with 1,500 nodes. For every value of a we report 
the average of 100 random realizations. First, we analyze the goodness of the upper bound 
found here for the average path length. In Fig. |6](a) we illustrate the variation of (/) with 
the connection radius for two different values of a, namely a = 5 and a = 10. In the same 
plot we illustrate the values of the upper bound obtained with the expression flTT)) . where 
it can be seen that the upper bound is very close to the average shortest path obtained 
for these RRGs. Particularly, for large values of r the observed values are almost identical 
to those of the upper bound. As observed in Fig. [6](b) the average path length not only 
changes with the variation of the connection radius but also with the rectangle side length. 
This is already expected from the eq. ffTSj) where it can be seen that as a —)■ cxo the average 
path length also grows to infinity for fixed r (see observation at the end of the previous 
paragraph). 

These results agree with our intuition that as we elongate the rectangle there are nodes 
which are farther apart from each other and as a result the average path length of the whole 
graph increases. For a better analysis of this relation we plotted the results of the average 
path length for 100 random realizations of the previously mentioned RRGs versus a, with 
m = 11, 250 edges, in Fig. |6](c). We also plot here the upper bound (ITT)) . It can be seen that 
there is an almost linear increase of (/) for values of 1 < a < 15 after which the dependence 
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(a) (b) (c) 

Figure 6. (a) Plot of the average path length for RRGs with n = 1,500 nodes with different 

connection radii for a = 5 (black squares) and for a = 10 (red squares). The solid squares represent 
the average of (1) for 100 random realizations (lines connecting the squares are used to guide the 
eye). The dotted lines represent the upper bounds (I17p for (1) in RRGs for the corresponding values 
of a. (b) Change of the average path length versus radius for networks with n = 1,500 nodes for 
different values of a: black squares (a = 1), red circles (a = 5), green triangle (o = 10), blue 
rhombus (a = 30). (c) Illustration of the upper bound (solid line) for the average path length as 
a function of a for a network with the same size as before and connection radii that guarantees 
m = 11, 250 edges for a given value of a. 


is very flat. In this region we have that a —?■ oo, which corresponds to a good approximation 
of a one-dimensional RGG. For a = 1 it is known that the average path length depends only 
on the inverse of the radius, (Ig) = 0 (l/r). The actual radius used for the plot in Fig. [7] 
(c) is r = 0.0578, which gives an estimate of the average path length of 17.3, which is not 
too far from the observed value in the plot for a = 1. In the case of a = 30 we are in the 
presence of a very elongated rectangle, which is very similar to a one-dimensional RGG. A 
crude estimate of the average path length in this case would be (le) = n/ (k), which in the 
current case will give (Ig) ~ 100, which is relatively close to the observed value of (/) 75 

for a = 30. 
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E. Clustering coefficient 


One of the most important network parameters is the local clustering coefficient, usually 
mown as the Watts-Strogatz clustering coefficient, of a node i. This parameter is defined as 

3: 

2h 


C,= 


(19) 


ki{ki - 1 )’ 

where t* is the number of triangles involving the node i and ki is the degree of the node 

i. Taking the mean of these values as i varies among the nodes in T, one gets the average 

clustering coefficient of the network: (C) = ^ 

The average clustering coefficient of a RGG has been obtained by Dali and Ghristensen 
„ log n -\- a 

17IJ for r = -—- when n —)■ cxo and a ^ oo, where a G M is a constant for a given 


riTT 


number of nodes (see Section B. Gonnectivity): 


{C,) = 


1 — Hfi (1) d even 
|-ffrf(l/2) d odd, 


( 20 ) 


where d is the dimension of the hypercube in which the nodes are embedded and 


d/2 


3\*+2 




where T (i) is the Gamma function. Thus, for d = 2 , {C 2 ) = 1 — 


3\/3 

47r 


( 21 ) 


0.5865 and for 


d=l, {Cl) = 3/4 = 0.75. 

Here, however, we are interested in an expression that accounts for the variations of the 
clustering coefficient with both the connection radius and the rectangle side length. Our 


strategy is similar to the one used in 


17| . That is, let i and j be two connected nodes 


m 


a RRG, which are separated at a Euclidean distance 5 from each other. Let us draw two 
circles of radius r centered respectively at i and j. Let 6 < r such that the two nodes are 
connected. Then, because 6 < 2r the two circles overlap. Because i and j are connected, 
any point in the area formed by the overlap of the two circles will form a triangle with the 
nodes i and j . In addition, any node inside the two circles which is not in the overlapping 
area forms a path of length two with the nodes i and j . Thus if we quantify the ratio of 
the overlapping area to the total area of the circle we account for the ratio of the number of 
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triangles to open triads in which the nodes i and j take place, i.e., the clustering coefficient. 
This ratio is given by 


2 r^ arccos 


{C) = 


— 1 - -(5V4r2 - 52 
2 r / 2 


7rr^ 


( 22 ) 


At this point we only need an estimation of the length 6 between two connected nodes 
in a RRG. We use here a simple approach based on the following intuition. Let us start by 
considering n nodes in a square in such a way that they form a regular square lattice. Then, 
5 is proportional to the length of the side of the rectangle a divided by the number of circles 
along this side. As we have a square, the number of points along the side of length a is the 
same as that for the other side. Consequently, 5 ~ If we elongate the rectangle to 

a ^ 00 , which resembles a straight line, we will have that the separation between the two 
points is just the length of the straight line divided by the number of nodes, 6 ~ an~^. For 
a general rectangle the separation between two points in a line along the edge side of length 
a is given by (j~ an~'^, where 7 ~ a (a + b)~^. Notice that when a = 1 (6 = 1) we have 
5 ~ and when a —?■ cxd (6 —)■ 0) we have 5 ~ an~^. We remind the reader that here we 
consider only the case b = a~^. 

Plugging S ~ an~'^, with 7 ~ a (a+ 6 )”^ into fl22l) we obtain some surprising results. 
In Fig. El^a) we illustrate the dependence of (C) with a for different radii based on the 
equation fl2^ with the estimated value of 6 given before. Notably, the clustering coefficient 
is predicted to change non-monotonically with the rectangle side length. Instead, for small 
values of a the clustering coefficient is predicted to increase to a maximum value and after 
it the clustering decreases linearly. In addition, according to this model, as the connection 
radius increases the clustering coefficient is expected to increase for the same value of a. In 
closing, the largest value of the clustering is expected for certain specihc value of a and a 
relatively large connection radius. 

In order to verify these hndings we compute the average clustering coefficient of RRGs 
with n = 1, 500 and the same connection radii used in the simulations with the analytical 
formulas. The results of the variation of the clustering with the rectangle side length are 
illustrated in Fig. d^b). As can be seen the clustering increases up to a maximum, whose lo¬ 
cation depends on the connection radius, and then decays with the increase of the elongation 
of the rectangle. We have not been able to capture the dependence of 6 with the radius in 
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our previous reasoning, but we very well captured the behavior of the clustering of having a 
non-monotonic change with a. Also, these experiments show that the increase of the connec¬ 
tion radius increases the average clustering coefficient as predicted by our analytical results. 
As can be seen in Fig. 0 (b) for a = 1 and small radius the average clustering coefficient is 


{C) 

to 


0.61, which is very close to the expected value for the 2-dimensional RGG according 


ITj. When a = 30and the radius is relatively large, the average clustering coefficient is 
~ 0.75, which coincides with the exact value expected for the one-dimensional RGG 


(G) 
according to 


m Gonsequently, the RRG generalizes the values of the clustering coefficient 


of both, the one- and two-dimensional RGG, for a = 1 and a —)■ oo, respectively. In addition, 
it provides a series of intermediate values of the clustering coefficient for intermediate values 
of the side length of the rectangle. 




(a) (b) 

Figure 7. Illustration of the dependence of the clustering coefficient with the rectangle side length for 
different connection radii: r = 0.05 (black), r = 0.1 (red), r = 0.15 (green) r = 0.20 (blue). In plot 
(a) we show the analytical results and in plot (b) the observed ones. Both results are obtained for 
RRG with 1,500 nodes and for the observed ones we averaged the results of 100 random realizations. 

We now further explore the relation between the radius r and the clustering for RRGs 
with different side lengths. We consider graphs with n = 1,500 nodes and a = 1,5,10,30. 
As the radius increases the graph is becoming more and more dense, which makes that the 
clustering coefficient is characterized by an abrupt increase at the beginning of the plot and 
then a linear increase until the value of {C) = 1 is reached for the complete graph (see Fig. 

E]). 
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Figure 8. (color online) Variation of the average path length with the radius for RRGs with 1,500 
nodes and a = 1, 5,10,30. Every point is the average of 100 random realizations. 

6 CONCLUSIONS AND FUTURE OUTLOOK 

We have introduced here a generalization of the RGG in which we embed the points into 
a unit rectangle instead of on a unit square. We consider a rectangle with sides of lengths 
a and h = 1/a, such that as when a = 1 we have the particular case of the classical random 
geometric graph embedded in a unit square. Also, when a —?■ cxd we have a very elongated 
rectangle which resembles a one-dimensional RGG. We have provided computational and 
analytical evidence that reaffirm the fact that the topological properties of the RRG differ 
significantly from those of the RGG. In this respect we have obtained analytical expressions 
or bounds for the average degree, degree distribution, connectivity, average path length and 
the clustering coefficient of RRGs. In general, these properties depend on the connection 
radius as well as on the rectangle side length. Most of the dependencies found here for these 
properties in terms of the rectangle side length are monotonic. The only exception is the 
clustering coefficient. This index first increases up to a critical value of a which depends 
on the connection radius, and then decays linearly with the increased elongation of the 
rectangle. 

The introduction of the RRGs opens new possibilities for studying spatially embedded 
random graphs. For instance, the study of dynamical processes taking place on the nodes and 
edges of these graphs is of great interest to explore how the shape constraints influence the 
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dynamics on the RRGs. On the other hand, the analysis of rectangular proximity graphs, 
such as the rectangular Gabriel graphs and random rectangular neighborhood graphs is 
also interesting for many of the practical applications of these graphs as mentioned in the 
Introduction. The generalization of the RRG model to higher dimensions is also of both 
theoretical and practical interest. In closing, the current work is expected to open new 
horizons for the study of random spatial graphs and its applications in physics and beyond. 
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